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ABSTRACT 

We have conducted a sensitive water maser search with the Australia Telescope Com- 
pact Array towards 267 1.2-mm dust clumps presented in the literature. We combine 
our new observations with previous water maser observations to extend our sample to 
294 1.2-mm dust clumps, towards which we detect 165 distinct water maser sites to- 
wards 128 1.2-mm dust clumps. Within the fields of our observations, we additionally 
find four water masers with no apparent associated 1.2-mm dust continuum emission. 
Our overall detection rate of 44 per cent appears to vary as a function of Galactic 
longitude. We find that there is an excellent correspondence between the locations of 
the detected water masers with the peak of the target 1.2-mm dust clump sources. As 
expected from previous similar studies, the water masers are chiefly detected towards 
the bigger, brighter and more massive 1.2-mm dust clumps. 

We find further evidence to suggest that the water masers tend to increase in flux 
density (and therefore luminosity), as well as velocity range, as the sources evolve. We 
also show that the current sample of water maser sources suffer less from evolutionary 
biases than previous targeted searches. 

We have compared the locations of the water masers with dust clumps which have 
a previously determined association with 6.7-GHz methanol masers and 8-GHz radio 
continuum. We flnd that a higher fraction of 1.2-mm dust clump sources in our sample 
are only associated with water masers (41) than only associated with methanol masers 
(13). This suggests that water masers can be present at an even earlier evolutionary 
stage than 6.7-GHz methanol masers. Comparison of the water maser detection rates 
associated with different combinations of methanol maser and radio continuum, as 
well as those with neither tracer, shows that the highest detection rate is towards 
those sources which also exhibit methanol maser emission. 

We have tested a previously hypothesised model for water maser presence towards 
1.2-mm dust clumps. We find water masers towards a high proportion of the clumps 
that the model predicts would have associated water masers but also find a number 
of water masers towards sources with a low calculated probability. We propose that 
this is likely an artifact of the poorly determined distances to the sources. We suggest 
refinements and future work which will further constrain the nature of the driving 
sources associated with water masers. 
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1 INTRODUCTION 

To date there have been relatively few large-scale, unbi- 
ased surveys for water masers. Instead, searches have fo- 
cused on high-mass star formation regions, typically tar- 
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geting re gions se l ected on the basis of IRA S colou rs (e.g. 
Churchw ell et al.l Il99(]| : ICodella etal] Il995: Suna da et all 
2007111 or othe r maser species (e.g iBatchelor et alTri980l : 



Caswell et al]|l983bl : ISeuther et afl 120021 '). Due to this, the 



known sample of water masers lacks the completeness that 
has been achieved f or both 6.7-G Hz methanol masers (e.g. 
ICaswell et al., ,201Q : iGreen et al] |201Q) and main-line OH 
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masers (e.g. lCaswel3ll998l ). This, coupled with the fact that 
the positions of the majority of southern hemisphere wa- 
ter masers are not known to sub-arcsecond accuracy, means 
that relatively little is known about the nature of sources 
that are able to produce luminous water masers. 

Large-scale, unbiased surveys for water masers are hin- 
dered by the large amounts of time required, and the sus- 
ceptibility of their observation to poor weather conditions. 
Further complications arise from the intrinsic variability of 
the sources which requires high positional precision in ei- 
ther the initial observations, or rapid follow up observa- 
tions. Some syst ematic searches for wat e r masers have been 
attempted (e.g. [Matthews et al.l Il985l: iBreen et al.l 12007 



that produce luminous water masers and the evolutionary 
phase they trace. Here we perform the necess ary observa- 
tions to test and refine the model developed in lBreen et al.l 



Hill et al.l l|2005l ) present a catalogue of 404 1.2-mm 



^- - . - , J 12003; 

ICaswell fc Breed I2OI0I : ICaswell. Breen fc EllingsenI I2OI1I ') 
but generally cover only very small regions of the Galactic 
plane. However, the observations of a relatively shallow, but 
large-scale single-dish survey for water masers have recently 
been completed, coverin g 100 degrees in l ongitude of the 
Southern Galactic Plane ( Walsh et al.ll20()8l ) and the survey 
results are expected in the coming year. 

Due to the difficulties associated with complete 
searches, attempts have been made to model water maser 
presence towards high-mass star formation regions, with 
the aim of gaining a more complete sample of water maser 
sources. These attempts have used particular observable 
properties of the high-mass st ar for mation regions to make 
these predictions. P alla ot al] (|l99ll ) modelled the probabil- 
ity of water maser presence on the b asis of IRAS f ar-inf rared 
colours. This model was tested by IPalla et all (|l993l ') and 
resulted in a detection rate of around 5 per cent. This dis- 
appointing result is most likely due to the poor spatial res- 
olution of the 5 observations (30 arcseconds at 100 ^m) 
which results in significant source confusion in crowded high- 
mass star formation regions; demonstrating the need for 
high-resolution complementary data in the formulation of 
such models. 

A systematic search for water masers of almost half 
a square degree o f the G 333.2-0.6 giant molecular cloud 
l|Breen et al.l 120071 ) revealed a strong relationship between 
water maser presence and the bigger, brighter, more 
massive and dense 1 .2-mm dust clumps (identified by 
iMookeriea et al.] (|2004l )), within the giant m olecular cloud. 
From these observations. iBreen et all l|2007l ) produced what 
appears to be a reliable statistical model for predicting which 
1.2-mm dust clumps will and will not have an associated 
water maser. While the likelihood of water maser presence 
increases with increasing values of all dust clump proper- 
ties tested, the simplest model with the greatest predictive 
properties was found to only include the dust clump radius. 
This model predicts that all 1.2-mm dust clumps of radius 
1.25 pc (or higher) have a probability of 0.5 (or higher) of 
having an associated water maser. 

While the model presented in iBreen et al.l (|2007l ) is 
promising, the water maser survey observations on which 
this model was derived, were of comparatively low sensi- 
tivity (detection limit of ~5 Jy or more in some regions). 
This, combined with the fact that the sample size was small, 
means that its validity cannot be determined before testing 
and refining it on a much larger sample. A reliable model 
relating the properties of 1.2-mm dust clumps to the pres- 
ence/absence of water masers represents a unique tool for 
investigating questions such as the mass range of the stars 



dust clumps, a large but manageable sample, perfect for 
a targeted wate r maser search to test the current model. 
iHill et all (2005) targeted their 1.2-mm dust continuum ob- 
servations towards 131 regions that were suspected of un- 
dergoing massive star formation, using the presence of pre- 
viously identified methanol masers and/or UCHll regions to 
select their targets. The observations were carried out using 
the SIMBA instrument on SEST (the 15 m Swedish ESQ 
Submillimetre telescope) and detected emission directly as- 
sociated with all but 20 of th e methano l mase rs and nine of 
the UCHll regions targeted. iHiU et all (|2005l ) also made a 
large number of serendipitous detections within the target 
fields. 

Here we present new Australia Telescope Co mpact Ar- 
ray (A TCA) observations towards 267 of the iHill et ahl 
(|2005h dust clumps and supplemen t our observations 
with data taken from the literature ([Breen et al.l l2010bl : 
ICaswell fc Breenll2010l : iForster fc Caswelllll989h . Altogether 
we present water maser data toward s 294 dust clump. Us- 
ing the model of iBreen et al.l l|2007l ) we find that few of 
these dust clumps have high calculated probabilities of wa- 
ter maser presence (only 58 have probabilities greater than 
0.01 1. 

iBreen et al.l ()2010ah present 12.2-GHz methanol maser 
observations targeting a number of these dust sources. This 
search for 12.2-GHz methanol masers, c onstitutes a near - 
complete sample within the sample of iHill et al] (|2005l ) 
sources, since no strong 12.2-GHz emission is expected in 
regions lacking 6.7-GHz methano l maser emission . Since a 
number of the sources presented in lHill et al.l (|2005l ) are con- 
tai ned within the regions searched for OH maser emission 
bv lCaswelil (|l998l ). a large subset of the 1.2-mm dust clumps 
have been searched for OH masers. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

The feasibility of this study was tested with observations us- 
ing the University of Tasmania's 26 m Mount Plea sant radio 
telescope. All 404 1.2-mm dust clumps given in et al.l 
(2005) were targeted in these observations during 2007. An 
estimated 150 water maser sources were detected within the 
regions. An effort was made to locate the origin of the emis- 
sions by observing four-point grids around each of the tar- 
gets and fitting a 2D-Gaussian to the relative amplitudes of 
the emission detected in each. Given the immense success 
in detecting water maser sources, together with the impos- 
sible task of confidently disentangling nearby sources and 
associating them with the 1.2-mm dust clumps (especially 
given the limited pointing accuracy of the instrument), we 
proposed to observe the sources with the ATCA. Using the 
ATCA has several advantages, foremost of which are the 
ability to obtain accurate positions and gain much greater 
sensitivity. 

Successful ATCA observations were carried out in three 
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Table 1. ATCA observations: epochs, array configurations and 
typical synthesised beam sizes. 



Array 


Epoch 


Beam size 


H75 


2007 October 9, 10 


28x24 


H214 


2007 December 7 


13x8 


H214 


2008 July 12, 13 


13x8 


6B 


2008 August 15, 16, 17 


1.7x0.5 



different array configurations over four epochs; these are 
hsted in Table [T] F or 137 of t he 40 4 1.2-mm continuum 
sources observed bv iHill et all (|2005l ). observations at 22- 
GHz are not presented. The reasons for the omission of 
these 137 sources were chiefly due to their: proximity to a 
dechnation of zero degrees (where observations can be par- 
ticularly troubling at the ATCA); location in the northern 
hemisphere; lack of reported dust characteristics; or because 
the weather conditions they were taken in were too poor to 
result in reliable data. Observations of an additional 27 1.2- 
mm dust clumps were not necessary since they are associated 
with water masers reported to good positional accuracy in 
the literature. 

In total, we successfully searched 267 1.2-mm dust 
clumps for water maser presence with the ATCA. For all 
observations the field of view of the telescope is almost 5 
arcmin, and the full width half maximum of the primary 
beam is 2.3 arcmin. Observations were carried out in groups 
of sources in the same neighbourhood so as to decrease over- 
heads and allowing reference pointing observations to be car- 
ried out efficiently. Pointings were carried out on a strong 
continuum source, located within ~20 degrees of the target 
sources, approximately every hour. Using this method, the 
pointing is usually accurate to 5 arcsec, compared to about 
20 arcsec when no pointing observations are carried out. 

Targets were observed in a series of cuts, interspersed 
with phase calibration observations. At all epochs of obser- 
vations only a single linearly polarised signal was recorded. 
The general observation strategy differed depending on the 
array and weather conditions, but in general, nearby sources 
were grouped together and bracketed by phase calibration 
observations carried out on a compact continuum source, lo- 
cated not more than 10 degrees from the targets. The time 
spent observing the targets and the time between phase cali- 
bration was tailored to most efficiently observe the targets in 
the allocated time on the ATCA. The average strategy was 
to observe each target field for 10 minutes, over five cuts, 
with phase calibration (onsource for 1.5 minutes) carried 
out every 10 minutes. In general, these observations have 
allowed 5-a" detections to be made at the 150 - 200 mjy 
level in the central section of the beam. For some sources, 
integration times were adjusted depending on the scheduled 
time allotment, and due to this the sensitivity limit may be 
either slightly lower or slightly higher (these are listed in 
column 10 of Table El. 

On each observing day, observations of a bandpass cal- 
ibrator were completed. Two different sources were used; 
PKS B1253-055 and PKS B1921-293, both are strong at 
22 GHz and are usefully observable at high elevations near 
the beginning and end of the time that the Galactic plane is 
visible, respectively. Primary fiux calibration is with respect 



to observations of PKS B1934-638, which were also carried 
out daily. 

Details of the the observations carried out in each of 
the array configurations are given below. Descriptions of the 
correlator configurations (and therefore spectral resolution 
and velocity coverage) and individual observation strategies 
are also given in this section. 



2.1.1 H75 

The first series of observations were carried out in the most 
compact ATCA configuration, H75. Primarily, observations 
of sources located near declination zero were intended for 
this array configuration, but we were able to observe some 
additional, early rising sources in a Director's time allocation 
prior to our scheduled time on 2007 October 9. 

For these observations the correlator was configured to 
record 1024 channels across a bandwidth of 16 MHz for a 
single polarisation. This correlator configuration allowed for 
a channel spacing of 0.21 km s~^, corresponding to a velocity 
resolution of 0.25 km s^^. The velocity coverage of these 
observations was just over 200 km s~^. 

In these observations, target sources were observed for 
10 minutes each, over a series of five cuts (in general) spread 
over several hours. Observations of a phase calibrator source 
were carried out every 8-14 minutes. 



2.1.2 H214 

Water maser observations were carried out in the H214 hy- 
brid configuration over two epochs. The first epoch was 
during 2007 December 7 and included only one source 
(G 213.705-12.597) which was observed as part of a mul- 
tifrequency study of the masers associated with nearby star 
formation region Mon R2. The adopted correlator configura- 
tion for these observations were very different for this source, 
using a bandwidth of 4 MHz and 1024 channels. In this con- 
figuration the velocity coverage is limited to ~50 km s~^, 
but the spectral resolution is increased to 0.06 km s~^. 

The main series of observations carried out in this array 
configuration were completed in 2008 July. For these obser- 
vations the correlator was configured to record 512 chan- 
nels over a bandwidth of 32 MHz for a single polarisation. 
We made this change in correlator configuration, as from 
our first series of observations, it was clear that for some 
sources we were unable to observe the full extent of the ve- 
locity range of the emission. In these observations we made a 
sensible compromise, forgoing additional spectral resolution 
for twice the bandwidth. The resultant velocity coverage is 
over 400 km s~^ and the spectral resolution is 1.0 km s~^. 
The average strategy for these observations was to observe 
each source for 12-14 minutes over 6-7 cuts of 2 minutes 
each spaced over several hours. Phase calibrator observa- 
tions were completed approximately every 13 minutes. 



2.1.3 6B 

The final series of observations was carried out in a 6B ar- 
ray in 2008 August. The correlator was configured as for the 
observations carried out in the H214 array, giving a velocity 
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coverage of more than 400 km s~ , and a spectral resolu- 
tion of 1.0 km s~^. Similar to the other observations, the 
observing strategy for observations in this array saw sources 
observed as a series of 2 minute cuts spaced over several 
hours. On average, sources were observed seven times, giv- 
ing onsource integration times of 14 minutes. Observations 
of a phase calibrator were conducted approximately every 
seven minutes. 



2.2 Data reduction 

All ATC A data were reduced using t he miriad software 
package (|Sault. Teuben fc Wrighj[l995l 'l. applying the stan- 
dard techniques for ATCA spectral line observations. Image 
cubes of the entire primary beam and usable velocity range 
were produced for each source. The flux densities of sources 
that were located away from the centre of the primary beam 
have been corrected to account for beam attenuation using 
the MIRIAD task 'linmos'. Both vector and scalar average 
spectra of the calibrated u^^data were inspected for each of 
the targets and cross checked with the emission identified 
in each of the image cubes. Spectra for each of the detected 
sources were produced by integrating the emission in the im- 
age cubes. The typical resultant rms noise in each spectrum 
was 25 - 40 mjy for sources located near the centre of the 
beam. 

The typical synthesised beam sizes experienced in each 
of the array configurations are listed in Table [T] In both of 
the hybrid configurations (II75 and H214), data from the 
6 km antenna (antenna 6) was not included in the data 
reduction. The adopted water maser rest frequency was 
22235.07985 MHz. 



3 RESULTS 



iHill et ahl (|2005l ) presented a catalogue of 404 1.2-mm dust 
clumps. Towards 294 of these clumps, we present water 
maser data and compare our detections with the derived 
probabili ty of water maser presence (using the model pre- 
sented in iBreen et al.1 (|2007l )). Most of the water maser data 
presented are new observations (132 of 165 water masers) , 
however, where a. p propr iate, data fromlBreen et al.l (l2010bl) . 
ICaswell fc Breed |2010') and iForster fc Caswelll l| 19891 ') has 
been used to minimise repetition of observations. 

Table [5] shows the targeted 1.2-mm dust clumps, along 
with their calculated probability of water maser presence, 
followed by the water maser data. We find a total of 165 wa- 
ter maser sources towards 128 of the 1.2-mm dust clumps. 
Spectra of each of the 132 water masers detected in our cur- 
rent observations are presented in Fig. [T] and are essentially 
ordered as the 1.2-mm dust clumps presented in iHill et al.l 
l|2005h (in order of increasing right ascension) unless there 
was a need for nearby sources to be vertically aligned. For 
the majority of sources, a velocity range of 200 km s~^ is 
shown, however, this has been reduced to 100 km s~^ for 
sources observed in the II75 configuration of the ATCA 
(which were carried out with a smaller bandwidth) and for 
one source, G 213.705-12.597, a velocity range of 40 km s~^ 
is shown (this source was observed with 4 MHz bandwidth) . 
Confusion from strong nearby sources are marked on the in- 
dividual spectra, except where they are present as features 



of negative intensities. In some spectra, where confusing fea- 
tures from multiple nearby sources are present, we use an 
'X' to indicate that those features are not associated with 
the presented source. A consequence of the relatively coarse 
spectral resolution of our observations, combined with the 
often strong and relatively narrow water maser features, is 
that 'ringing' is present in a number of spectra. 

A number of the sources that we detect have 
been previously presented in the literature (e.g 



V 

Johnston et al. 19721: Caswell et al.l 1974 : Kaufmann et aU 



19761: ICenzel fc Downed Il977l: ^Batche lor et all 119801: 



Braz fc Scaliselll982l: iBraz fc Epch tcin 1981; ICasweU etall 
19891 : iHofner fc Churchwelll Il996l . and references therein). 
However, the majority of these observations were conducted 
20 or more years ago, with relatively poor positional 
accuracy. When combined with the often extreme variable 
nature of water maser sources, the task of identifying our 
sources with those in the literature is therefore fraught with 
difficulties. Since the only certain way to match up sources 
with those in the literature is by position, we only compare 
our detections with observations of comparable positional 
accuracy. 

Several of the sources that we detect ha ve been ob- 
served in interferometri c obse r vations (e.g. iBreen et al.] 
2010bl: ICaswell fc BreenI I2OI0I : ICaswell fc PhillipsI l200a 
Forster fc CaswelJ Il989iv allowing us to compare our po- 
sitions and therefore assessing our positional uncertainties. 
In general, our positional agreement with previous obser- 
vations is excellent, someti mes better than 0.5 a rcsec, but 
consistently within 2 arcsec. iBreen et all l|2010bD present a 
more extensive assessment of the positional uncertainties of 
their water maser observations carried out with the ATCA 
and conclude that perhaps the largest contributor to the un- 
certainty was the tendency for the water maser spots to be 
spread out and the difference between measuring a slightly 
different feature at a different epoch. Here we similarly find 
that this is the case, and therefore expect the positional 
uncertainties of the water masers that we present to be ac- 
curate to 2 arcsec. 

Another facet of this assessment allowed us to deter- 
mine the effect that the synthesised beam has on the posi- 
tional accuracy of the sources. As can be seen in Table [TJ 
the beam sizes resulting from the hybrid configurations are 
much larger than that of the more extended 6B array. Com- 
parisons with other measurements shows that even when the 
beam is large, accurate positions can be derived, since it is 
rarely the case that two sources show emission at the same 
velocities and are spatially separated by less than a synthe- 
sised beam. The much smaller beam of the 6B array has 
allowed for some water maser sites to be broken up into in- 
dividual components, although it is likely the case that very 
close companions are intimately associated with the same 
exciting source. 

Notes on some individual sources are presented in Sec- 
tion 13.11 We find that there is an excellent correspondence 
between the locations of the detected water masers and the 
1.2-mm dust sources and this is discussed further in Sec- 
tion 14.31 In addition to the water masers that we detect 
towards 1.2-mm dust clumps, we detect four water maser 
sources that seem not to be associated with any dust con- 
tinuum emission and these are discussed in Section [3.21 



Table 2. Target 1.2-mm dust clumps, probability of water maser associatio n, followed by a description of the water maser observations including any detections. Columns 1-3 give the 
1.2-mm du st clump name foll owed by the right ascension and declination l|Hill et al.ll2005l ) and column 4 shows the calculated probability of water maser presence (using the model 

Column 5 gives information about the ATCA array configuration used in the water maser observations followed by the epoch. In so me cases, data is 
which is marked in column 5 with a 'BCEP' followed by the year of the observations. This format is also followed for sources taken from lCaswell &: BreenI 



Breen et al. 



derived by 

taken from lBreen et al _ _ 

1I2OIC1I ) fCB) and iForster'^r^aawelll lll989l ) (FC). Columns 6-13 give details of the water masers detected towards the dust clumps (which in some cases is more than one source), as well 
as 5-(T detection limits where no source is detected, specifically; column 6: water maser Galactic coordinates; column 7; water maser right ascension; column 8: water maser declination; 
column 9: separation between water maser and 1.2-mm dust clump peak (arcsec); column 10: water maser peak flux density (Jy) or if preceded by '<' the 5-(7 detection limit (mjy); 
column 11; velocity of the water maser peak emission (km s^^); column 12: velocity range of either the detected emission, or the observed velocity range if no detection; and column 
13: integrated flux density of the water maser emission (Jy km s~^). Column 14 presents notes related to the sources, which in most cases are further discussed in Section |3. II (note the 
following abbreviations: BS(RS) - indicates that the source is dominated by a blue(red)-shifted feature.) 



Dust clump Water maser 

Name RA Dec prob. array Name RA Dec Sep. 

{l,b) (J2000) (J2000) and (Z,6) (J2000) (J2000) (") 

(degrees) (h m s) (° ' ") epoch (degrees) (h m s) (° ' ") 



Sp 

(Jy) 



Vp 



Vr 



Int notes 



G213 


61-12.6 


06 


07 


47 


9 


-06 


22 


57 


3 


7E-4 


H214dec07 


G213 


705-12.597 


06 


07 


47 


85 


-06 


22 


56 


6 


1 


92 


12 


4,13 


78 


G269 


45-1.47 


09 


03 


13 


5 


-48 


55 


22 




0.15 


H214jul08 


G269 


461-1.471 


09 


03 


15 


01 


-48 


55 


31 


4 


18 


0.6 


51 


49,52 


1.6 
























H214jul08 


G269 


457-1.467 


09 


03 


14 


87 


-48 


55 


12 


6 


24 


0.8 


66 


56,69 


1.7 


G269 


15-1.13 


09 


03 


32 


3 


-48 


28 


00 


1 


2E-4 


H214jul08 


G269 


153-1.128 


09 


03 


33 


52 


-48 


28 


03 





12 


26 


7 


6,16 


40 


G270 


254-0.84 


09 


16 


41 


4 


-47 


55 


46 


4 


lE-5 


H214jul08 


G270 


254-f0.835 


09 


16 


41 


13 


-47 


56 


11 


3 


25 


98 


9 


4,17 


199 
























H214jul08 


G270 


257-1-0.834 


09 


16 


41 


62 


-47 


56 


17 


4 


31 


60 


8 


-2,48 


227 


G284 


271-0.391 


10 


23 


47 





-57 


48 


38 


8 


5E-5 


H214jul08 
























<155 




-200,150 




G284 


307-0.376 


10 


24 


04 





-57 


49 


02 


1 


5E-3 


H214jul08 
























<155 




-200,150 




G284 


338-0.417 


10 


24 


06 





-57 


52 


06 


5 


4E-5 


H214jul08 
























<400 




-200,150 




G284 


35-0.42 


10 


24 


10 





-57 


52 


39 


1 


4E-4 


H214jul08 


G284 


351-0.418 


10 


24 


10 


92 


-57 


52 


33 


9 


9 


233 


7 


-32,61 


847 


G284 


345-0.404 


10 


24 


12 





-57 


51 


42 


1 


5E-4 


H214jul08 
























<280 




-200,150 




G284 


341-0.389 


10 


24 


14 





-57 


50 


46 


5 


8E-3 


H214jul08 
























<215 




-200,150 




G284 


328-0.365 


10 


24 


15 





-57 


49 


10 


4 


9E-5 


H214jul08 


G284 


329-0.365 


10 


24 


15 


27 


-57 


49 


09 


4 


2 


0.8 


1 


-5,1 


1.7 


G284 


384-0.441 


10 


24 


18 





-57 


54 


46 


1 


4E-4 


H214jul08 
























<185 




-200,150 




G284 


344-0.366 


10 


24 


21 





-57 


49 


42 


1 


5E-4 


H214jul08 
























<156 




-200,150 




G284 


352-0.353 


10 


24 


27 





-57 


49 


18 


1 


9E-3 


H214jul08 
























<200 




-200,150 




G287 


37-1-0.65 


10 


48 


05 


2 


-58 


26 


40 


3 


lE-4 


H214jul08 
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Table 2. - continued 



Dust clump Water maser 
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Table 2. - continued 
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Figure 1. Spectra of the 22-GHz water masers detected towards 1.2-mm dust clumps. 
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3.1 Individual sources 

Here we draw attention to sources that we were unable to 
adequately describe in Table [2] Close neighbouring water 
masers are discussed, as well as sources with surprising ve- 
locities and intriguing histories. For some sources, associa- 
tions with other maser species are also discussed. 

G 213.705-12.597. This water maser is located in a 
nearby star formation region, Mon R2. Associated with this 
star formation region are many maser species, including 6.7- 
and 12.2-GHz methanol, main-line OH, and some excited 
OH maser transitions. Of particular interest is the excited 
OH maser at the 4765-MHz transition; these sources are very 
rare (only about 20 known in the Galaxy) and this source is 
by far the strongest known at this transition. Furthermore, 
the emission from this 4765-MHz OH maser undergoes flar- 
ing ac tivi ty that has not been seen in any other maser of this 
type (jSm its 2003; Smits, Cohen & Hutawarakorn 1998). 

G 284. 351-0. 41 8. This source is also presented in 

iBreen et al.] ll2010bh as it is coincident with an OH maser 
l|Cas well"l998'). although is named G 284.350-0.418, due to 
a slight difference in mea sured position affecti ng the rounded 
Galactic coordinates. In iBreen et al.l l|2010bl ). spectra from 
two epochs are shown to display the typical level of vari- 
ability seen in the water maser sources. When the 2003 and 
2004 spectra are compared to our current observations, it 
can be seen that the source has once again varied by a mod- 
erate level and has increased in peak flux density by about 
a factor of 2. 

G 291.271-0.719, G 291.274-0.709 and G 291. 284- 
0. 716. These three sources are also listed in iBreen et al.l 
l|2010b[ ) since they all fall within the f ield of a target OH 
maser source. ICaswell fc PhiUipsI i 20081 ) also observed this 
group of sources, and remark that the strongest source, 
G 291.284-0.716, is a member of a distinct class of sources 
that are dominated by strong blue shifted outflows. In fact, 
maser emission has never been detected at the systemic ve- 
locity of the source (to a limit of 0.3 Jy in these observa- 
tions), which is around 100 km s~^ from the detected emis- 
sion. The location of G 291.284-0.716 is very close to the 
edge of the 1.2-mm dust clump emission, unlike the major- 
ity of sources. 

G 291.579-0.431 and G 291.579-0.434. G 291.579- 
0.431 shows emission spread over approximately 2 arcsec and 
in Table [2] we present the median of the positions measured 
for the individual features. This source is clearly separated 
from G 291.579-0.434, which is located more than 10 arcsec 
away. 



G 301.136-0.226. As noted in iBreen etHI (|2010bh . 
emission from this source is significantly spread. In the ob- 
servations carried out in 2008 July, we observe individual 
maser components spread over almost 5 arcsec. An OH 
mase r is located wi thin the spread of the water maser emis- 
sion l|Caswelllll99i ). 

G 305.362+0.150. During the observations in 2008 Au- 
gust we detected a moderately strong water maser exhibiting 



one ma in spectral feature. In 2004 (presented in lBreen et al.l 
(|2010bl )) an additional, broader feature of more than 100 Jy 
was also detected and has disappeared intervening interval. 

G 0. 546-0. 851b. This water maser falls within a cluster 
of three distinct water maser sites and has a velocity range 
that exceeds 170 km s~^. The emission extends almost sym- 
metrically about the peak water maser emission. 

G 5.513-0.255. Like G 291.284-0.716, this water 
masers lies unusually close to the edge o f the 1.2-mm dus t 
clump emission that was observed by iHiU et all l|2005h . 
No other similarities between the sources are evident. This 
source is very weak and shows emission close to the likely 
systemic velocity of the region, while G 291.284-0.716 is very 
strong and is a dominant blue-shifted source. 

G 8.727-0.395. This water maser is clearly projected 
against 1.2-mm dust continu um emission on inspection of 
the corresponding dust map of lHill et al.l l|2005l ), although no 
dust clump is reported at this position. It is likely that this 
source was omitted from the catalogue of iHill et al.l (|2005l ) 



G 10.473+0.027 and G 10. 480+0. 034. The observa- 
tions pr esented in Table [2] have been taken from Bre en et al.l 
(|2010bD . Both of the sources were detected in observations 
carried out in 2003, but G 10.4804-0.034 had fallen below 
the detection limits in 2004 (<0.2 Jy). 

G 11. 903-0. 142. The observations presented in Ta- 
ble[2]have been tak en from the observations carried out in 
iBreen et al] (|2010bl ). and similarly to G 10.480-1-0.034, this 
water maser was detected in 2003, but by 2004 the emission 
had fallen below our detection limits. Since the emission de- 
tected in 2003 had a peak flux density of 0.3 Jy, this is not 
surprising. 

G 12.209-0.101, G 12.429-0.048, G12.720-0.218 and 
G 13.874+0.281. These four well separated water maser 
sources all exhibit water maser emission at only negative 
velocities. From the velocities of nearby methanol masers, 
systemic velocities in this region of the sky would be ex- 
pected to be between 20 - 60 km sT^ . Therefore these sources 
are all potential candidates for the class of sources that are 
dominated by blue-shifted emission. 



3.2 Water masers with no apparent associated 
dust continuum emission 

In addition to the water maser sources presented in Table [2] 
and Fig. [TJ we detect several other sources which are of spe- 
cial interest. While an exhaustive attempt to search for wa- 
ter maser emission beyond the boundaries of the 1.2-mm 
dust clumps was not made, here we list the additional water 
maser sources that were identifled. 

Interestingly, we flnd four water maser sources that ap- 
pear not to be associated with any 1.2-mm dust cl ump emis- 
sion a nd these are presented in Table[3]and Fig.[2l lHill et al.l 
(|2005h detected 1.2-mm dust clump emission in the fields of 
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all of the target 6.7-GHz methanol masers and UCHll re- 
gions, but for 20 of the methanol masers and 9 of the UCHll 
regions there is no spatially coincident dust continuum emis- 
sion. Investigation by Hill et al. (2005, ) showed that these 
sources appeared to be no different to those where 1.2-mm 
dust continuum was detected. Since the presence of 1.2-mm 
dust clump emission indicates the presence of cold, deeply 
embedde d sources, t wo of the possible explanations put for- 
ward by iHill et alj (|2005l ) are that these sources are more 
evolved and hence no longer in the cold core phase, or, per- 
haps they are less massive. 

We have investigated possible mid-infrared counter- 
parts for these four sources by comparing their locations 
with products from the GLIMPSE survey. All are pro- 
jected against regions of extended infrared emission, and two 
sources (G19. 612-0.120 and G 23.455-0.201) are coincident 
to within 1 arcsec of a GLIMPSE point source. Interestingly, 



the two sources with associated point sources are also asso- 
ciated with methanol masers. This suggests that these two 
sources may be located at the far distance and consequently 
were not detected in 1.2-mm dust continuum due to sen- 
sitivity limitations. The nature of the other two sources is 
more confusing and certainly warrants future investigation. 

Unlike 6.7-GHz methanol masers, water masers have 
been detected towards both evolved and low-mass stars. 
Their presence in locations devoid of 1.2-mm dust contin- 
uum emission is therefore much more easily accounted for 
than methanol masers. Comments on each of these four wa- 
ter masers are given below. 

G 305.248+0.195. This water maser is separated from 
the boundary of the nearest dust clump by more than an 
arcminute. 

This source falls beyond the latitude range completely 
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searched for the prese nce of water maser emission by 
ICaswell fc Breed l|2010l ) which was focused on a region 
bounded by a fongitude range of 305.0 - 306.26 ° and ± 
0.15° latitude. A total of 23 water maser sources were de- 
tected in this portion of the Galactic plane. Analysis of the 
masers found in this area, along with another region (longi- 
tude range of 311.0 - 312.18° and ± 0.15° latitude) showed a 
surprising preponderance of sources with highly blu e-shifted 
emission. Investigation bv lCaswell fc BreerJ l|2010h led to a 
suggestion that these water masers represent a population 
of sources that are tracing a very early evolutionary stage, 
perhaps preceding the stage where 6.7-GHz methanol masers 
are seen. 

Interesting, G 305.248-1-0.195 only shows emission at a 
blue-shifted velocity of -95 km s~^ (nearby methanol masers 
detected in the Methanol Multibeam Survey exhibit peak 
velocities between -40 and -30 km s~^), and therefore adds 
to the already high number of blue-shifted sources in this 
region of the Galactic plane. However, the interpretation 
that such sources may be v ery young directly contradicts one 
of the suggestions given bv lHill et al.l (|2005l ) as to why some 
sources are devoid of 1.2-mm dust continuum emission - that 
they are more evolved. Alternatively, it is possible that this 
water maser is associated with a low-mass object, although a 
high velocity feature > 50 km s~^ from the systemic velocity 
of the region would be unusual for a low-mass star. 

G 19.612-0.120. This water maser is coincident with 
one of the methanol masers from IWalsh et al.] l|l998l ) that 
was targeted bv lHill et al] (j2005l ). but is also offset from the 
boundary of t he ne arest dust clump by about an arcminute. 
I Green et al.l (|2010l ) present evidence suggesting that this 
methanol maser is located at the far distance, which would 
limit the sensitivity of the 1.2-mm dust continuum observa- 
tions, but not to the extent that no detection of dust con- 
tinuum emission would be expected. 

G 23.397-0.219 and G 23.455-0.201. Both of these wa- 
ter masers fall within the one dust continuum emission 
map, but neither are directly associated with any detectable 
dust continuum emission. The first source is apparently soli- 
tary (i.e. without any associated UCHll region or 6.7-GHz 
methanol maser emission), but the second sour ce appears to 
be coi ncident with an UCHll region reported bv I Walsh et al.l 
l|l998h . 



4 DISCUSSION 

4.1 Detection statistics 

Altogether, we find water maser emission towards 128 of the 
294 1.2-mm dust clumps, equating to a detection rate of 44 
per cent. Investigation of the detection rate of water masers 
towards 1.2-mm dust clumps shows some variation with 
Galactic longitude. The water maser detection rate towards 
dust clumps located between longitudes of 213 and 333 de- 
grees is 55 per cent (51 of 93), compared with 38 per cent 
(77 of 201) for the longitude range to 38 degrees. Although 
not covering exactly the same longitude ranges, it is curi- 
ous to note that the detection rates of 12.2-GHz methanol 
masers in the longitude range 270 to 305 degree s were signif- 
icantly lower than in other parts of the Galaxy (|Breen et al.l 



l201lh . If it were the relative distances that were affecting 
our detection rates then, if anything, we would expect to 
have a lower detection rate in the 213 to 333 degree longi- 
tude region as we would expect relatively more dust clump 
detections in these regions. Since our maser observations 
are of a high sensitivity we would expect a minimal detec - 
tion rate dependence on distanc e. Caswell fc Breer] (|2010l ) 
and ICaswell. Breen fc EUingsenl (|2011 ) completed system- 
atic searches for water masers within three regions of the 
Galactic plane. They show that there is a large difference in 
the water maser population densities within these regions, 
suggesting that local factors have an effect on the detection 
statistics. It is therefore possible that there are broader lo- 
calised effects contributing to our detection statistics here. 
This warrants further investigation, but is beyond the scope 
of this paper. 



4.2 Characteristics of the detected water masers 

4.2.1 Velocity ranges and flux densities of the complete 
sample 

The average velocity range of the 165 water masers is 26 
km s~^ and they have a median velocity range of 13 km s~^. 
In comparison, the ave rage velocity range of the 379 wa- 
ter masers detected by iBreen et al.1 l|2010bl ) is 29 km s"^ 
and they have a median velocity range of 15 km s~^. While 
these numbers are quite similar, it is intriguing that veloc- 
ity ranges of the sourc es that are p o tentiall y biased towards 
more evolved sources (|Breen et al.l (|2010bl ) chiefly targeted 
their observations towards OH masers) are, in general, larger 
than those detect ed towards the sample of dust clumps. 
iBreen et al.l (|201ll ) presents evidence that suggests that both 
the 6.7- and 12.2-GHz methanol maser sources show an in- 
crease in velocity range as the sources evolve. Although com- 
plicated by the presence of high- velocity features, we have 
carried out a similar analysis on the water masers in Sec- 
tion |3321 

Further evidence that these water masers are less biased 
towards more evolved sources lies in th e average and median 
flux density of the detected sources. iBreen et al.1 l|2010bD 
present evidence that the water masers increase in flux den- 
sity as the sources evolve. The average flux density of the 
water masers detected towards the dust clump sample is 
39 Jy and they have a median of 4 Jy which are both lower 
than the average (57 Jy) and median ( 5 Jy) fl ux density 
of the sources presented in iBreen et al.l (|2010bl ) . This sug- 
gests that in the dust clump sample, we have detected water 
masers over a broader range of evolutionary stages. 



4.2.2 Association with other maser species and 
continuum 



iHiU et al.1 l|2005l ) determined associations between their 1.2- 
mm dust clump sources and the 6.7-GHz methanol masers 
and UCHii regions (detected at 8 GHz) towards which their 
observations were targeted. We have made an investigation 
of the number of d ust clumps , with and without associated 
water masers, that lHill et al.l l|2005l ) have designated as be- 
ing: 'mm-only' cores (i.e. sources without associated radio 
continuum emission or 6.7-GHz methanol masers); associ- 
ated with 6.7-GHz methanol masers; associated with both 
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Table 3. Four water maser sources detected in the target fields that are not coincident with any detected 1.2-mm dust clump emission. 
Columns 1-7 give the: Water maser source name; water maser right ascension; declination; peak flux density (Jy); velocity of the water 
maser peak emission (km s""*^); velocity range (km s~^); and the integrated flux density of the water maser emission (Jy km s""*^), 
respectively. 



Name 


RA 


Dec 


Sp 


Vp 


Vr 


Int 


(l,b) 


(J2000) 


(J2000) 


(Jy) 


(km s-^) 


(km s^^) 


(Jy km s^^) 


(degrees) 


(h m s) 


(° ' ") 










G 305.248+0.195 


13 11 34.84 


-62 35 11.0 


7 


-93 


-95,-92 


9 


G 19.612-0.120 


18 27 13.55 


-11 53 14.7 


1.2 


58 


57,60 


1.6 


G 23.397-0.219 


18 34 42.32 


-08 34 43.3 


7 


102 


97,106 


14 


G 23.455-0.201 


18 34 44.89 


-08 31 07.4 


3.1 


82 


57,84 
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Figure 2. Water maser sources with no associated 1.2-mm dust clump emission, detected within the target flclds. 



6.7-GHz methanol masers and radio continuum; and finally, 
those associated with only radio continuum. In all of our 
statistical analysis, 293 of the 294 targets for water masers 
have been included (omitting a single source which is miss- 
ing values for some dust clump properties). 

Table |4] shows a summary of the categories of dust 
clumps where we detect the water masers, as well as the 
water maser detection rate in each category. We find an 
overall water maser detection rate of 44 per cent towards the 
targeted dust clumps, and this number falls to 23 per cent 
for the 'mm-only' dust clump sources. The highest water 
maser detection rates are towards those dust clump sources 
that exhibit methanol maser emission, closely followed by 
those dust clumps showing only radio continuum emission. 
It can be seen that there are many more dust clump sources 
showing only water masers (41), than those showing only 



methanol masers (13). While some fraction of these 1.2-mm 
sources that show only methanol masers are likely to have 
associated water masers that were below our detection limits 
at the time of observation, it seems unlikely that this could 
be the case for all of them (by considering the statistics of 
detectability for the two epochs of water maser observations 
presented in iBreen et al . (2010b)). However, these sources 
will provid e interesting ta rgets for further water maser ob- 
servations. iHill et al] l|2005i ) suggests that their 'mm-only' 
cores may represent an earlier evolutionary phase than those 
sources showing methanol maser emission. If this is the case, 
these numbers certainly imply that water masers can be 
present even earlier than methanol masers. A supporting 
argument for water masers having a longer lifetime than 
that of methanol masers is that water molecules are both 
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Table 4. Water maser dete ction rates in each dust clump cate- 
gory (from IhIU et alj l|200d )): 'mm-only' (those sources showing 
no associated methanol maser or radio continuum emission); as- 
sociated methanol masers (meth); associated with both methanol 
masers and radio continuum associations (meth-|-cont); and asso- 
ciated continuum (cont). These numbers have been presented in 
two categories; those dust clumps with no detectable water maser 
emission, and those where we have detected water masers. The 
fourth column shows the water maser detection rate (per cent) in 
each of the categories. 





No water 


Water 


Detection rate 
(%) 


'mm-only' 


137 


41 


23 


Meth 


13 


48 


79 


Meth-|-cont 


5 


25 


83 


Cont 


10 


14 


58 


Total 


165 


128 


44 



more easily created and more robust than methanol, which 
is consistent with our interpretation. 

An alternative explanation for the higher number of 
dust sources showing only water masers (than only methanol 
masers) is that a number of these additional water masers 
are associated with lower mass objects. Comparing the dust 
clump properties in each of the association categories of- 
fers no evidence for there being any systematic differences 
between the groups which may result if a large number of 
these sources were lower-mass objects. 

We have investigated the water maser characteristics 
in each of the dust clump association categories. Table [5] 
presents the average and median values of both the water 
maser peak flux density and the velocity r ange. Like the 
methanol masers presented bv lBreen et al.l l|201lh . we find 
that there is a general gradient which trends from lower 
values of both flux density and velocity range for the (prob- 
ably) younger sources to higher values for the more evolved 
sources. This trend progresses through the categories as fol- 
lows: those associated with 'mm-only' sources, to those with 
methanol masers, and then to those with methanol maser s 
and radio continuum (|Breen et al.ll2010al : [HiU et al.ll2005h . 
The exception to this trend is the last stage, where only 
radio continuum is seen. 

If the numbers presented in Table[5]are representative of 
the larger population of water masers, then they imply that, 
while the water masers increase in flux density and velocity 
range as the sources evolve, towards the end of their lifetime 
this trend 'turns over' and they steadily decrease in both 
flux density and velocity range, presumably decreasing until 
the maser emission ceases. This is contrary to that seen in 
the case of methanol masers (|Breen et al.ll201ll ). which seem 
not to experience any downwards trend, implying that they 
switch off much more abruptly. 

Comparing the luminosities of the wa ter m asers (calcu- 
lated using the near distances presented in lHill e t al. (2005)) 
with the velocity ranges of the water maser sources, sim- 
ilarly shows an increase in both luminosity and velocity 
range through the different association categories from least 
to most evolved. However, there is significant overlap be- 



Table 6. Slopes and intercepts of the line of best fit of log inte- 
grated maser luminosity versus log of the H2 number densities of 
the 1.2-mm dust clumps. Values for 6.7- and 12.2-GHz methanol 
masers are listed for comparison with the water maser values for 
comparison. Errors are presented in square brackets and follow 
the estimates of each value. 



Maser Slope Intercept Reference 
species 

6.7-GHz -0.72[0.25] 5.53(1.17] Breen et al. (2011") 

12.2-GHz -0.56[o.28] 4.16(1.23] Breen et al. C2010a 'l 

water -0.43(0.07] 3.86(0.33] current 



tween the values within each association category. A high 
level of scatter is expected when considering water masers, 
since a number of sources show high velocity features and it 
is not clear that these are restricted to a certain evolution- 
ary stage. Also, the sources that are associated with only 
8 GIfz radio continuum are scattered throughout the range 
of values that are associated with the sources in other cat- 
egories, but none have large values of luminosity, and very 
few have large velocity ranges. These findings are consis- 
tent with those conclusions drawn from the average values 
that are presented in Table [S] where peak flux densities are 
considered. 



4-2.3 Comparison between water maser luminosity and 
1.2-mm dust clump H2 number density 

iBreen et al.l l|2010al ) shows that the H2 number density of 
the 1.2-mm dust clum ps decreases as the sources evolve. 
Although, as noted by iBreen et al.l l|2010ah . the apparent 
change in density might be a con sequence o f the constant 
temperature assumption applied bv lHill et al.l l|2005l ) in their 
calculations and consequently the difference in densities 
might instead represent an increase in temperature. In either 
case, the changes to the physical properties of the source are 
consistent with evolution. 

In Fig. [3] we present the log of the integrated water 
maser luminosity versus the log II2 number densities of the 
associated 1.2-mm dust clumps. T his plot show a slope that 
is consistent with the findings of iBreen et al.l l|2010ah and 
those previously presented in this section - that the more 
luminous water masers are more evolved and that these are 
associated with the less dense 1.2-mm dust clumps and vice 
versa. Table |6] compares the computed slopes and intercepts 
of the line of best fit of the log luminosity 6.7- and 12.2-GIIz 
methanol masers, along with the water masers versus the H2 
number densities of the associated clumps. 

Inspection of the values listed in this table show firstly 
(as previously stated) that the 6.7-GHz methanol masers 
increase in flux density more rapidly as they evolve than 
the associated 12.2-GHz sources, and, as shown by the 
larger intercept, are generally stronger than the 12.2-GHz 
sources. Interestingly, the values associated with the water 
maser sources imply that the water masers also increase in 
flux density less rapidly than the 6.7-GHz methanol masers 
as they evolve and are also weaker in general at a given 
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Table 5. The average and median flux densities, and velocity ranges of all of the detected water masers. Water maser sources have been 
broken up into several categories, and are as for Table [J] with the addition of 'all sources' which shows the values for the entire sample. 



Water 


Average 


Median 


Average 


Median 


classification 


flux density 


flux density 


vel range 


vel range 




(Jy) 


(Jy) 


(km s-^) 


(km s-^) 


all sources 


39 


4.4 


26 


13 


'mm-only' 


2 


1.5 


19 


8 


Meth 


44 
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28 


19 


Meth+cont 


107 


18.5 
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1 


26 


17 



"8 





O o 


Oo 





















o o 




ooo o 







o 




o ° 


o ° ^ 








o 
o 



o o 



o 

O 00 

o o o 



o o 



~~1 I I I 

3.5 4.0 4.5 5.0 

Log H_2 number density 



5.5 



Figure 3. Log water maser integrated luminosity (Jy km s~^ 
kpc^) versus log of the H2 number density (cm~^) of the associ- 
ated 1.2-mm dust clump. The dashed line shows the line of best 
fit (which has a slope of -0.43[0.07] and an intercept of 3.86[0.33] 
and a correlation coefficient of 0.50). 



the observations presented here were completed only once, 
we cannot directly determine the level which maser variabil- 
ity has affected our results. Instead, we derive predictions 
about the completeness of these o bservations by comp aring 
this sample with the f indings of iBreen et al.l (|2010bh and 
ICaswell fc BreenI (I2OI0I). 



Breen et al.1 (|2010bl ) finds that 17 per cent of the 253 wa- 



ter maser sources that were observe d at two epochs we r e only 
detectable at one of these epochs. ICaswell &: BreenI (|2010l ) 
found similarly that 16 per cent of the 32 sources that they 
detected in their double epoch complete searches were de- 
tectable only once. If we simply apply the larger percentage 
(i.e. 17 per cent) to our current sample, it might be expected 
that ~34 additional water maser sources would be detected 
if a second epoch of observations were conducted towards 
the target 1.2-mm dust clumps. Since we detect 165 water 
masers towards 128 1.2-mm dust clumps, we would perhaps 
expect that these additional sources would be detected to- 
wards a further 26 1.2-mm dust clumps. 

However, the situation is a little more complex than 
this. It is not clear if the role of variability would result in 
many more dust clumps being recognised as being associ- 
ated with water masers, or alternatively that different water 
maser sources would be detected at different epochs towards 
the same dust clump. It is likely that both would be true, 
that is, that some number of water masers would be detected 
towards additional 1.2-mm dust clumps and that some ad- 
ditional/different water masers would be detected towards 
1.2-mm dust clumps where we have detected water masers. 



II2 number density (than the sample of 6.7-GHz methanol 
masers). However, it is difficult to determine which of the 
maser species are present at the earliest evolutionary stage of 
high-mass star formation and therefore which species shows 
stronger emission at a given evolutionary stage. 



4.2.4 Water maser variability 

Water masers are known to ex hibit extreme var iability over 
relatively short timescales (e.g.lFelli et al.l 2007). The water 
maser observations presented in .Breen et al.l ( 201 Obi ) were 
completed, in many cases, twice with the o bservation epochs 
separated by ~10 months. Additionally, ICaswell fc BreenI 
(|2010l ) carried out a double epoch search for water maser 
emission of two small regions of the Galactic plane. Since 



4.3 Location of the water masers in the 1.2-mm 
dust clumps 

Analysis of the locations of the water masers wi th respect 
to the 1.2-mm dust clump peak, as reported by iHill et al.l 
(|2005D . shows that there is excellent correspondence be- 
tween the ir locations. The 1.2-mm dust continuum maps of 
iHill et al] (2005) have 8 arcsec pixels, and the listed peak po- 
sitions are just the central position of the peak pixel. On av- 
erage, we find that the separation between the water maser 
positions and the associated dust clump peak is 13 arcsec 
and that the median separation is 8 arcsec. It is therefore 
likely that, in the majority of cases, the water masers are 
very closely associated with the peak 1.2-mm dust clump 
emission. 

Fig.|4]shows a histogram of the separations between the 
dust clump peak and the detected water masers. All water 



Water masers towards 1.2-mm dust clumps 25 




2 4 6 8 10 12 14 16 18 20 22 more 



Separalion between dust peak and water maser (arcsec) 

Figure 4. Histogram of the angular separation between tlie re- 
ported 1.2-mm dust clump peak and the position of tlie water 
maser sources. 



masers that were separated by more than ~30 arcsec from 
the dust clump peak were further investigated by comparing 
the reported dust clump peak positi o n with the dust maps 
presented in figure Al of iHill et ahl (|2005l ) and the water 
maser positions. We find in just over half of these cases (11 
of 19) there is an obvious error in the reported dust clump 
position which has lead to an apparently large separation 
from the water maser position. In the worst examples, the 
reported dust clump peak is located beyond the bou ndary 
of the respective clump. Even though [Hill et af] (|2005l ) pub- 
lish images of their 1.2-mm dust clumps, the resolution (and 
contrast) is insufficient to enable us to determine better posi- 
tions independently. Given the obvious problems with some 
number of the dust clump peak positions which erroneously 
result in some of the most extreme separations between the 
peak of the clump and the associated water maser, it is 
likely that the actual average and median separations be- 
tween the peak of the dust clumps and the associated water 
masers would be even smaller, perhaps comparable to the 
pixel size. 

We have compared the locations of the water masers 
with respect to the location of their associated dust clump 
peak and the peak flux density of the maser sources. It is 
evident that very few of the water masers that are truly 
significantly offset from the peak of the dust clumps have 
flux densities that are greater than a few Jy. The possible 
corollary, that most of the water masers that are located far 
from the dust clump peak are relatively weak, needs further 
investigation. 

4.4 Assessing the current model of water maser 
presence in 1.2-mm dust clumps 

iBreen et all (|2007l ') created a model for water maser pres- 
ence towards 1.2-mm dust clumps which was derived 
from observations of th e 1.2-mm dust continuum emission 
l|Mookeriea et al1l2004l l and water maser observations that 
they completed within the G333 giant molecular cloud. By 
determining which 1.2-mm dust clumps were associated with 
water masers a s well as those tha t were devoid of water 
maser emission, IBreen et al.1 ||2007D fitted a Binomial Gen- 
eralized linear model to the maser presence/absence data, 
using the 1.2-mm dust clump properties as predictors. They 



found that the simplest model with the greatest predictive 

properties was based only on dust clum p radius. 

We have used the model derived bv iBreen et al.1 (I2007D 
to calculate the probability of wate r maser d e tectio n for 
each of the dust clumps presented in iHiU et al.1 (|2005l ) that 
have dust clump radi us measurements. For a number of dust 
sources presented in iHill et al] l|2005h there is a near-far 
distance ambiguity. In these cases, we have assumed the 
near distance for c onsistency with analysis completed by 
iBreen et al.l (|2010al ). However, we note that in this instance 
this assumption might be particularly hazardous, since the 
predictive model relates only to 1.2-mm dust clump radius 
which is highly influenced by the adopted distance measure- 
ment. 

Using the model for water maser presence associated 
with 1.2-mm dust clumps, it was determined that only 58 
of the 294 dust clumps considered in this analysis, had a 
probability of 0.01 or greater of having an associated water 
maser. Probabilities for individual dust clumps are presented 
in Table[2] Comparison between the calculated probability of 
water maser and the actu a l dete ctions shows that the model 
described in iBreen et al.1 l|200'iD is promising. The average 
probability of water maser presence in the 128 1.2-mm dust 
clumps where water masers are detected is 0.10, and has a 
median value of 0.00081. In comparison, for the dust sources 
where we detect no water maser emission, the average prob- 
ability of water maser presence is 0.027 and has a median 
value of 0.00026. An alternative description of these statis- 
tics is that, of the 27 1.2-mm dust clumps which had a calcu- 
lated probability of > 0.1, 20 yielded water maser detections. 
If the considered probability is lowered to be >0.01 then the 
number of predicted clumps is increased to 58 and the num- 
ber of clumps with detections is 39. It is therefore clear that 
the model is picking up a high detection rate for the sources 
with high calculated probabilities, however, given that we 
detect water masers towards 128 dust clumps, it is also clear 
that the model is failing to return high probabilities of wa- 
ter maser presence in around two thirds of the clumps with 
detections. 

Given that the original model was derived from a sam- 
ple of water masers that were detected in a complete search 
with a relatively poor sensitivity limit (>5 Jy), the obvious 
property to investigate is the relative flux density of the de- 
tected water masers. It might be expected that there would 
be a trend whereby the water maser detections towards the 
1.2-mm dust clumps with high association probabilities have 
higher flux densities, given the nature of the data on which 
the model was derived. The average flux density of the water 
maser sources detected towards dust clumps with a proba- 
bility of 0.01 or higher (58 water maser sources) is 48 Jy and 
these sources have a median flux density of 6 Jy. In compar- 
ison, the average flux densities of the water masers detected 
towards clumps with a lower calculated probability is 33 Jy, 
with a median flux density of 3.5 Jy. It is possible, therefore, 
that the current model is biased towards predicting correctly 
the presence of the stronger sources. However, more likely, is 
that the sources with very low probabilities are dominated 
by far distance sources which accounts for the lower average 
flux density (which would also mean that their probabilities 
have been under estimated). 

Another obvious factor to consider is the high degree 
of dependence on distance measurements. Since the origi- 
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nal model was derived from sources located within the one 
GMC with a well constrained distance measurement, the 
poorly determined distances towards most star formation re- 
gions did not affect the generation of this model. The model 
for water maser presence associated with dust clumps uses 
only the radius of the dust clump in order to calculate the 
probability. As previousl y discussed , we h ave assumed the 
near distance reported bv lHill et al.l l|2005l ) for sources with 
distance ambiguities. It is clear that this is an unreliable as- 
sumption, especially since our search is of high sensitivity 
which reduces the likelihood of any significant detection de- 
pendence on distance. It would be expected that the split of 
near and far distance allocations should be approximately 
even (or maybe even favour the far distances due to the 
large volume). Inspection of the characteristics of some of 
the dust clump sources that have low calculated probabili- 
ties, yet have associated water masers, shows that if the far 
distance was used the calculated probability would be high 
(in many cases >0.3). The increase in radius that results in 
a change from near to far distance does not result in a high 
probability in water maser presence in all of the dust clump 
sample, especially those with no detected associated water 
masers. 

If our assertion, that the poorly constrained distance 
measurements are adversely affecting the predictive capa- 
bilities of the model, is correct, then the distribution of the 
clump properties should support this. If, for example, we in- 
spect the peak flux density of the 1.2-mm dust clumps with 
associated water masers and calculated probabilities of > 
0.01, they should not be any different from the sample of 
1.2-mm dust clumps with associated water masers and cal- 
culated probabilities of <0.01. But, they are in fact different. 
The average peak flux density of the clumps with calculated 
probabilities of >0.01 (and have associated water masers) is 
2.8 Jy beam~^, whereas the average of those clumps with 
lower than expected calculated probabilities (<0.01) is 1.4 
Jy beam~^. Again, the most likely explanation for this is 
that a number of the 1.2-mm dust clumps which have low 
calculated probabilities, but yet have an associated water 
maser, are located at the far distance. This means that the 
seeming downfall of the model to allocate high probabilities 
to a number of the 1.2-mm dust clumps where we detect 
water masers may be eradicated if accurate distance mea- 
surements were obtained. 

4.5 Fitting Binomial general linear models to the 
new water maser and 1.2-mm dust clump data 

While the current model for water maser presence in 1.2- 
mm dust clumps is promising, we ha ve repeated the analy- 
sis carried out by [Breen et all (|2007l ) with the much larger 
sample presented here. However, two important considera- 
tions are: firstly, the adoption of the near distance measure- 
ments for the sources with distance ambiguities; and sec- 
ondly, that the 1.2-mm dust clumps in this sample were not 
de tected in a systematic manner (unlike the observa tions 
of iMookeriea et al] (|2004l ') used in iBreen et al.1 (|2007h ). al- 
though given the high number of additional 1.2-mm dust 
sources found within the fields of the target sources this is 
probably not detrimental to the analysis. 

We have fitted a Binomial generalized linear model to 
the water maser presence/absence data, using the 1.2-mm 



Table 7. Analysis of deviance table for the single term models 
(using the 1.2-mm dust clump properties from lHill et al.l |2005)), 
showing the deviance and the AIC together with the associated 
likelihood ratio statistic and p-value for the test of the hypothesis 
that the stated single model provides no better fit than the null 
model consisting only of an intercept. 



Predictor 


Deviance 


AIC 


LRT 


p- value 


none 


399.84 


401.84 






integ 


341.81 


345.81 


58.02 


2.590e-14 


peak 


332.98 


336.98 


66.85 


2.923e-16 


fwhm 


354.23 


358.23 


45.61 


1.445e-ll 


dist 


398.65 


402.65 


1.19 


0.2761 


mass 


366.41 


370.41 


33.43 


7.391e-09 


radius 


379.44 


383.44 


20.40 


6.299e-06 


density 


399.66 


403.66 


0.18 


0.6697 



dust clump properties as predictors. Fo r a more deta i led de - 
scription of this analysis method, see iBreen et al. I (|2007l ). 
All but one of the targeted 1.2-mm dust clumps were in- 
cluded in this an alysis. The sourc e that has been excluded 
is G 6.60-0.08 as iHill et all (120051 ) suggest that the derived 
mass of this source is uncharacteristic of high-mass star for- 
mation regions. For each of the 1.2-mm dust clumps, all 
of the derived clump properties were tested: the integrated 
fiux density (Jy), peak flux density (Jy beam~^), source full 
width at half maximum (FWHM) (arcsec), distance (kpc), 
mass (Mq), radius (pc) and H2 number density (cm~^). In 
the following sections we consider p-values of less than 0.05 
to be statistically significant (i.e. the hypothesis that the 
single term model provides no better fit than the null model 
consisting only of an intercept is rejected when the p- value 
is less than 0.05). 

Firstly, all dust clump properties were tested to see if 
individually they could give an indication of the likelihood 
of associated water maser presence. This was done by fit- 
ting a single term addition Binomial model to each prop- 
erty (or 'predictor') and showed an increasing probability of 
the presence of water masers associated with increasing val- 
ues of dust clump integrated flux density, peak flux density, 
FWHM, mass and radius (shown in Table [7]). This result is 
simil ar to that foun d when investigating methanol masers 
( Breen et al.ll20l"0al ). and means that any one of these dust 
clump properties can give an indication of the likelihood of 
flnding an associated water maser. Box plots of each of the 
dust clump properties broken up into the categories of 'n' 
and 'y', referring to those clumps with no associated water 
maser, and those with an associated water maser, respec- 
tively, are presented in Fig. [S] These box plots show graph- 
ically the same information as the results of the single term 
additions; that lower values of the dust clump properties; 
integrated and peak flux density, FWHM, mass and radius, 
are associated with the 1.2-mm dust clump sources with no 
associated water maser emission, compared to those associ- 
ated with water maser emission. 

After determining the predictive capabilities of individ- 
ual dust clump properties, stepwise model selection was used 
to create the simplest model with the greatest predictive 
power. This method considers all of the dust clump prop- 
erties and tries to maximise the predictive properties while 
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Figure 5. Box plots of the 1.2-mm dust clump properties split 
in the categories of yes and no, referring to the presence of water 
maser emission. There is a statistically significant difference be- 
tween the two categories in the 1.2-mm dust clump integrated and 
peak flux density, mass, FWHM and radius. This information is 
the graphical display relating to the information in Table [7] Note 
that source G 23. 960-f 0.137 has been removed from these plots, 
since is is more than twice as dense as all of the other sources 
shown on the plot and therefore is an extreme outlier. 



Table 8. Summary table for the Binomial regression model, 
showing for each predictor the estimated coefficient and the stan- 
dardised z-value and p-value for the test of the hypothesis that 
ft=0. 



Predictor 


Estimate 


Std. Error 


z value 


p- value 


Intercept 


-1.563746 


0.271202 


-5.766 


8.12C-09 


peak 


0.994999 


0.254413 


3.911 


9.19e-05 


fwhm 


0.011744 


0.007242 


1.622 


0.105 



ous dataset where a model based on dust clump radius (pc) 
was found to be best, the current dataset uses dust clump 
peak flux density and FWHM (arcsec). As can be seen in 
Table [S] the dust clump peak flux density is the most influ- 
ential term, indicated by the large standard error relative to 
the coefficient of the dust clump FWHM together with the 
large p- value. 

The misclassification rates for this new model are some- 
what promising when the probability threshold is set to 
0.5, it correctly predicts half (63 of 127) of the dust clump 
sources with associated water masers, and it correctly pre- 
dicts 144 of the 166 dust clumps with no associated water 
masers. It is easy to understand why such a model would fail 
on a high number of the dust clumps with associated water 
masers, since we know that these sources have higher peak 
flux densities on average, but as can be seen in Fig. (5] there 
is a large overlap in the range of values of this property be- 
tween those clumps with associated water masers and those 
without. This is consistent with a number of dust clumps 
with associated water masers having apparently low peak 
flux densities, since they are located at the far distance. 
Therefore a peak flux density measurement that is scaled 
with reliable distance measurements (a pseudo luminosity) 
may offer a promising model. Likewise, a reliable model may 
be derived from radius measurements that have correctly ac- 
counted for distances. 



4-5.1 A new model based on 1.2-mm dust clump radius 

Since the model for water maser presence presented in 
iBreen et al.l (|2007h uses only the dust clump radius to pre- 
dict the probability of associated water maser presence, for 
comparison, a model using only dust clump radius has been 
produced from this new data. The resultant estimated re- 
gression relation is 



trying to minimise the number of terms (and therefore dust 
clump properties). The resultant model contains only two 
dust clump properties; peak flux density and FWHM. The 
estimated regression relation is 

log— ^ = -1.564 + 0.995a;peak + 0.012a;FWHM, 

where speak is the peak flux density, sfwhm is the source 
FWHM, and pi is the probability of flnding a water maser 
towards the i*** 1.2-mm dust clump. The regression summary 
of this model is shown in Table [S] In contrast to the previ- 



log-^^ = -0.981 + 1.45a;radius, 
1-Pi 

where a^-adius is the radius of the 1.2-mm dust clump, and pi 
is the probability of finding a water maser towards the j*'' 
1.2-mm dust clump. The regression summary of this model 
is shown in Table (9] 

Since only one dust clump property is present in the re- 
gression relation, it is easy to determine the physical impli- 
cations of the model. For example, if the probability of water 
maser presence is set to 0.5, the corresponding dust clump 
radius is 0.68 pc. This means that the model is saying that 
1.2-mm dust clumps with a radius of 0.68 pc (or higher) have 
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Table 9. Summary table for the Binomial regression model using 
only dust clump radius, showing the estimated coefficient and the 
standardised z-value and p-value for the test of the hypothesis 
that /3i=0. 



Predictor Estimate Std. Error z value p-value 

Intercept -0.9814 0.2101 -4.670 3.01e-06 

radius 1.4460 0.3544 4.081 4.49e-05 



a probability of 0.5 (or higher) of having an asso ciated water 
maser. In comparison, the model presented bv lBreen et al.l 
predicts that 1.2-mm dust clumps with radii greater 
than 1.25 pc for the same probability threshold of 0.5 will 
have associated water masers. The model produced with the 
current data, predicts that water masers are associated with 
much smaller dust clumps (about half the size; which is per- 
haps not surprising given that these observations are of much 
higher sensitivity), albeit relatively poorly. The misclassifi- 
cation rates from this new model are poor, particularly in 
determining which clumps will have associated water masers 
(which it gets wrong two-thirds of the time). 

The poorly constrained distance measurements are 
therefore having a large effect on our ability to model the 
water maser presence within the 1.2-mm dust clumps. How- 
ever, the higher accuracy of the models in predicting the 
sources without associated water masers correctly, may in- 
dicate that there is a population of sources that are signif- 
icantly far from the properties of the clumps that exhibit 
water masers that not even a change from near to far dis- 
tance would boost them to a high probability of water maser 
presence. Therefore, from the current data we have been un- 
able to derive an adequate model for water maser presence 
from dust clump radius. 



5 CONCLUSION 

We find water maser emission towards 128 of the 294 1.2-mm 
dust clumps searched. In total, we detect 165 distinct water 
maser sites and most are either new, or achieve accurate 
positions for the first time. 

There is an excellent correspondence between the posi- 
tions of the water masers and the peaks of the 1.2-mm dust 
clumps in the majority of sources. In addition to the wa- 
ter masers that we detect towards our target 1.2-mm dust 
clumps, we detect four sources towards regions apparently 
devoid of dust continuum emission (from comparison be- 
tween the water maser locations and the 1.2-mm continuum 
maps of iHill et al. ( 2005. )). We suggest that two of these 
sources may be located at the far distance which has re- 
sulted in lower sensitivity 1.2-mm continuum observations 
tow ards these s ources . 

iHill et al.l l|2005h allocated their 404 1.2-mm dust 
clumps in their catalogue the following categories: associ- 
ated with a 6.7-GHz methanol maser; associated with both 
a methanol maser and 8-GHz radio continuum; associated 
with radio continuum only; and millimetre only, i.e. not asso- 
ciated with either methanol masers or radio continuum. We 



have assessed our detection rates of water masers towards 
the dust clumps in each of these categories. The highest wa- 
ter maser detection rates are towards dust clumps which are 
associated with methanol masers (both with and without ra- 
dio continuum), and the lowest detection rate is towards the 
millimetre only sources (although at 23 per cent, this detec- 
tion rate is still quite high). We additionally find that there 
are more dust clumps that are only associated with water 
masers (41) than are only associated with methanol masers 
(13). This suggests that water masers can be present at an 
even earlier evolutionary stage than methanol masers. 

Comparison between the 1.2-mm dust clump properties 
with and without asso ciated water maser emission, shows 
(similarly to results of iBreen et al.l (|2007l )) that the water 
masers are associated with the bigger, more massive sources 
with higher peak and integrated flux densities. There is a 
trend whereby the more luminous water masers are associ- 
ated with 1.2-mm dust clumps with lower H2 number densi- 
ties than the less luminous water masers. This trend is also 
seen in the case of both 6 .7- and 12.2-GHz methanol masers 
(jBreen et al.ll2010al . [2oTll) . 

Like the 6.7- and 12.2-GHz methanol masers presented 
bv lBreen et"aLl l|201lh we flnd that there is evidence for both 
the luminosity and the velocity range of the water masers 
to increase as the sources evolve. This implies that the gas 
volume conducive to the maser emission also increases with 
evolution. The water maser sources show evidence for this 
trend to 'turn over' towards the end of their lifetime, pre- 
sumably showing a decline in both luminosity and velocity 
range until the emission ceases. 

We have used these water maser observations to test 
the model for wate r maser pres e nce t owards 1.2-mm dust 
clumps presented in iBreen" et al.l l|2007D . We find that there 
is a large number of water maser detections towards dust 
sources for which the computed probability of water maser 
presence in greater than 0.01, with a detection rate of 67 
per cent towards these sources. However, since the number 
of clumps where we detect water maser emission (128) is 
higher than the number of clumps for which the probability 
of water maser presence is greater than 0.01 (58) it is clear 
the the model needs some refinement. The inadequacy of the 
current model may be (at least partially) attributed to the 
adoption of the near kinematic distance for sources where 
distance ambiguities exist, since an assignment of the far 
distance would result in a larger probability in sources where 
the water masers are detected. 

We have attempted to create a new model for water 
maser presence towards 1.2-mm dust clumps, but find that 
we are severely limited by distance uncertainties. Our anal- 
ysis shows that the prospects of creating a reliable model 
for water maser presence within 1.2-mm dust clumps is high 
when reliable distances can be assigned to the so urces. Com- 
paring the success of our model with that of iPalla et al.l 
(jl991 ) shows that our model has much greater success in 
predicting water maser detectability. 

A crude evolutionary implication of our derived model 
(in conjunc t ion w ith the original model we present by 
iBreen et al.l (|2007l )). is that dust clumps with radii equal to 
0.97±0.29 pc (calculated using the difference between the 
radii threshold implied from each model) have a 50 per cent 
chance of forming one or more sources that are able to ex- 
cite water maser emission. Although a simplistic view, this 
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further implies that the lifetime of the dust clump is approx- 
imately twice that of water masers. 
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